Narrow entanglement beats 
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We study how the entanglement between two atoms can be created or modified even when they 
do not interact but when each of them interacts dispersively, i.e., weak and far from the resonance 
with a single mode of the field. Considering that regime we apply a method which makes use of a 
small nonlinear deformation of the usual SU (2) algebra in order to obtain the effective Hamiltonian 
describing correctly the dynamics for any initial states. In particular we study two cases: In the 
first one we consider each atom initially in a pure state and in the second case we assume that they 
start in a Werner state. We find that both atoms can reach, periodically, maximum entanglement 
if each of them starts in any eigenstate of a x , independent of the initial Fock state of the mode. 
Thus we find that a dispersive vacuum can generate entanglement between two two-level atoms. In 
the second case and when the field mode is initially in a coherent or thermal state, we find that in 
the high energy limit, in general, there is no entanglement between the two atoms however at well 
defined moments the initial entanglement is as suddenly recovered as removed. This time behavior 
looks like narrow beats separated by the so called entanglement dead valleys. 

PACS numbers: 03.67.-a, 03.65.-w 



I. INTRODUCTION 



II. THE HAMILTONIAN MODEL 



In 1935 E. Schrodinger introduced into the quantum 
world the entanglement concept by means of his com- 
munication addressing the gedanken experiment known 
as Schrodingers cat In the same year A. Einstein, 
B. Podolsky, and N. Rosen argued the incompleteness of 
quantum mechanics by describing the reality they sensed 
Q. Later, in 1964, J. S. Bell reported that a special 
no local operator must satisfy in average an inequality 
which can be violated only by some non separable states 
0. Today the non locality effect or entanglement is con- 
sidered as a resort for manipulation of quantum informa- 
tion which does not have a classical counterpart [J, Q. 
Thus, during the last two decades a major research ef- 
fort has been conducted in the emerging field of quantum 
information theory Q based on the renewed non locality 
effect. With this motivation, there has been a lot of in- 
terest in understanding and quantifying entanglement of 
pure and mixed states @, H, M, El • The entanglement of 
two systems can arise through direct interaction between 
them as well as through the coupling of the systems with 
a common quantum bus in the form of an auxiliary sys- 
tem or environment [ill, 03, 03 0, EE Ell • 

In this work we investigate whether and how the en- 
tanglement between two atoms can be generated or mod- 
ified when they interact dispersively with the same single 
mode of the electromagnetic field. We use the method of 
the small nonlinear deformation of the usual 5/7(2) alge- 
bra [13, [H, E3l > m order to obtain an effective Hamilto- 
nian describing the dynamics of the system. 



We consider two noninteracting two-level atoms, la- 
belled by sub or supraindexes a and b, each one coupled 
dispersively with a single mode characterized by the fre- 
quency lo. The unitary dynamics in the whole tensorial 
product Hilbert space, H = H <g> ® H mo d c , is driven 
by the Hamiltonian (h = 1) in the rotating wave approx- 
imation, 



H = 



» 



(b) 



+g a (a^b + aL a V) + g b (a^b + ^V), (1) 
where b and W are the annihilation and creation single 



mode operators, = |1 
is the z— component of the effective angular spin-half 
operator whose eigenstates are {|0)j, |l)j}, for j = a, b. 

Taking into account that the excitation number op- 
erator N = 
tion, [H,N[ 
H — loN + 



cr U) = \0)jj(l\, and a 



(f) 



with A„ 



= (o'i a ' ) + cri b - ) )/2 + b^b is a constant of mo- 
|=0, the ([T]) Hamiltonian can be written by 
Hint where 



r O) 



-gaia^b + cr^tf) 



g b (afb + a W b% (2) 



oj a — uj and A b — ui b — u>. 
We have assumed dispersive interactions between each 
atom and the common single mode. In other words, those 
couplings are weak and far from the resonance, so we can 
define the small parameters: 

< 7j=F= < 1, j = a,b, (3) 

*1 V( n )T 
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where (ti)t is the average photon number. Making use of 
the small rotation method [13, EH to obtain the effective 
Hamiltonian which approximately describes the interac- 
tion process, we can eliminate the two terms which do not 
represent the resonance interaction but represent rapid 
oscillations in the rotating frame. That can be achieved 
by applying to the Hamilnonian ^ a small unitary trans- 
formation: 



R 



tf)+e b (cr)°>b- 



Considering terms up to first order in e a and ej, of the 
Cambell-Baker-Hausdorf expansion, RHi nt B) , we obtain 
the following effective Hamiltonian: 



H, 



A„ 



v. 



(4) 
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The first two terms in the above effective Hamilto- 
nian represents the free evolution of the non-resonant 
atoms with renormalized transition frequencies. The 
third term is the so-called, dynamical Stark shifts, which 
describes an additional intensity dependent detuning of 
non-resonant atoms from the mode frequency. The last 
term represents an effective dipole-dipole interaction be- 
tween the non-resonant atoms which appears as a con- 
sequence of a collective nature of the interaction of non- 
resonant atoms with the quantized mode. We point out 
that this kind of effective interaction could not appear in 
the classical field and that the contribution of this term 
strongly depends on the internal resonance conditions of 
the non-resonant atoms. 

The effective Hamiltonian (0| can be diagonalized 
without difficulty but, for the sake of simplicity, we sup- 
pose that both atoms are identical in a way such that 
g = g a = g\y and A = A a = A&. Under those conditions 
the unitary evolution operator is given by 

u = e - l i A+92 ( 2fctb+1 )/ A i t |i) a |i) fca (i| fc <i| 

+e-"~< 2t / A \^ + )(4> + \ + e^'^I^XVn 

+el [A + ^(2 fc t b+1)/A ] t|0)a|0)ba(0|b(0|5 (5) 

where 1^) = (|0) a |l) 6 ± \l) a \0) b )/V2 are two Bell states. 

From Eqs. ([3]) and (O one realizes that the 
dressed states |l) a |l)b|n), |0) a |0)b|n), and |V' ± )|«) (n — 
0, 1, 2, . . . ) are stationary. Thus, by their form, indepen- 
dent of the initial mode state, the a — b bipartite sys- 
tem does not evolve when starting at one of the states 
|l) n |l)&, |0) a |0)b, or |i/' ± ), and hence preserves the initial 
entanglement. 



III. BIPARTITE ENTANGLEMENT 

The entanglement between two systems in a whole 
pure state is given by the entropy of any subsys- 



be evaluated by the concurrence 
cr„|^>*)|, where the asterisk de- 



tem. It can also 

notes complex conjugation of the probability ampli- 
tudes in the a z ® ^-representation, i.e., in the base: 
{|0)|0),|0)|1),|1)|0),|1)|1)} @,[Io|. The generalization of 
the concurrence to a mixed state p of two atoms is de- 
fined as the infimum of the average concurrence over all 
possible pure state ensemble decompositions of p, defined 
as convex combinations of pure states Si = {pi, IV'i)} de- 
composition, such that p — J2iPi\i ; i}(' l Pi\- I n this way, 
C(p) = inf Si ^iPiCdtpi}). Williams Wootters succeeded 
in deriving an analytic solution to this difficult minimiza- 
tion procedure in terms of the eigenvalues A^'s of the non- 
Hermitian operator pa y ®a y p* a y ®a y , where the asterisk 
again denotes complex conjugate of the elements of p in 
the a z ® <7 z -representation. The closed-form solution for 
the concurrence of a mixed state of two atoms is given 
by C{p) — max{0, a/AT — a/A^ — \f%i — a/AT}, where the 
Ai's are in the decreasing order [lOj. 

First we study the entanglement of formation between 
the two atoms when each of them is initially in a pure 
state. After that we consider the two atoms in a mixed 
states. In both cases we study how the entanglement 
evolves as a function of time and energy mode, under 
the weak and dispersive interaction ((4}. We consider 
the initial mode state to be a Fock state |n), a coher- 
ent state |a), and a thermal state px = Sn( n )T/(l — 
(n) T ) n+1 |n)(n|, where (n) T = l/(e"/ fesT - 1) is the av- 
erage photon number of the mode, T the absolute tem- 
perature and kg the Boltzmann constant (h = 1). 



A. Initial pure state 

Here we consider each atom to be initially in a pure 
state, i.e., qubit a being in \ip) a and qubit b in \ip)b- When 
the single mode is initially in a Fock state \n), the two 
atoms-system does not tangle with it at all time, and 
they evolve to the following pure state: 

|0) = {0\^){0\ip)^ A+ i( 2n+1 ^\0) a \0) b + \L) 



(i|^)(i|^) e - l [ A+ ^( 2 " +1 )] t |r i ,ii 



where \L) = Lo|0) a |l)b + Li|l) a |0)b is an unnormal- 
ized state in the subspace spanned by {|0) a |l)&, |l) a |0)b}. 
Here we have defined the functions: 

Lo = (0|^)(l|^)cos^- l (l|^)(0|^)sin^, 
La = mm cos £^_i<0|$<%) sin ^. 

Thus, the concurrence of the above pure state is read as 
follows 

C(\4>)) = 2|L ii - <0|V;)((%)<1|^)(%)| . (6) 

From Eq. ©, we can see that the concurrence: i) 
does not rely on the n photon number, ii) reaches the 
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maximum value 1 at t = 7rA/(2<? 2 ) and hence peri- 
odically, under the condition that both atoms start in 
any eigenstate of a x or even in a general state of type 
\9) = (|0) +e t6 \l))/V2 being 9 real. Therefore two atoms 
starting in any \9) state reach maximum entanglement 
when they interact dispersively even with a common vac- 



uum. In this form we can say that the vacuum can gener- 
ate entanglement between two two-level atoms even when 
they are far from the resonance. 

When the single mode is initially in a coherent state 
the reduced density operator of the two atoms- 
system becomes, 



p = |(l|^)(lM| 2 |l) a |l) ba (l| b (l| + \L)(L\ + |(0|^(0|^)| 2 |0) a |0) ba (0| b (0| 
+ (0|^(0|^)e^ 2 / A )*e-l"l 2 ( 1 - e2s! ' 2t/A )|0) a |0) b (i| 
+ (0|^)(V.|l)<0|^)^|l) e 24 ( A +^/ A )*e-l"l 2 ( 1 - e4ls2t/A )|0) Q |0) bQ <l| b <l| 
+ (,M0)^|0)e-^ 2 / A ) t e-l"l 2 ( 1 ^ 2l92t/A )|i) a <0| b (0| 
+ (^|l)Ml)e^ 2 / A )*e-l"l 2 ( 1 - e2l92t/A )|L) Q (l| b (l| 

+(V'|o)(i|^)(^|o>(i|^) e - 24 ( A +^ 2 / A )* e -i a i 2 ( 1 - e ^ ! ' 2t/A )|i) a |i) ba (o| b (o| 

+ (l|^(l|^) e - 4(A+s2/A)t e-l Q l 2 ( 1 ^ 2l92t/A )|l) Q |l) b (L|, (7) 



From this expression, Eq. ([7|), we see that only the 
last six non diagonal terms depend on the \a\ and they, 
in general, vanish for \a\ ^> 1. However, at times 
t = t k = nkA/g 2 (k — 1,2,...) those terms suddenly 
reappear and become independent of the intensity \a\. 
Since at times tk the state described by Eq. ||7J) is equal 
to the state when a = 0, that is, with an initial vacuum 
state, then at those times the concurrence is given by the 
Eq. ([5]) evaluated at tk which is zero. We can also see 
that when the atom a starts in the |0) state and the atom 
b begins in the |1) state or vice versa the concurrence of 
mixed state ([7]) is given by | sin(2<7 2 t/A)| reaching the 
maximum value 1 at t = irA/Ag 2 and repeating it peri- 



odically. On the other hand, in the high intensity regime, 
\a\ 1, the reduced density operator at any t ^ tk, is 
given by the first three diagonal terms of the Eq. J7]) and 
its concurrence is read as follows: 

C(p) = max{0,2(|L ii| - I W><0|v><lM<%>l)}- (8) 

In this regime and when each atom begins in any \9) 
state, the entanglement between them is always zero. 

When the initial mode state is a thermal state at abso- 
lute temperature T, the reduced density operator of the 
two atoms, at time t becomes: 



p = |(l|^)(l|^)| 2 |l) a |l) ba (l| b (l| + \L)(L\ + |(0|V)(0|^)| 2 |0) Q |0) ba (0| b (0| 

+ C ' (A+ ^ )t . 2g2t (W)(0|y)l0) g |0)»(£l + M1>M1>|L) <1| 6 <1|) 

l + (n) T (l-e^) 

+ e-( A +^)* {m{ i Ml)amL \ + M0>fo>|0>|L> a <0| 6 <0|) 

l + (n) T (l-e-^) 

(0|^)(0|y,)(^|l)(y,|l)e 2 ^+^) t 
l + (n) T (l-rf) 

(l|^)(l|^)(^|0)(^|0)e- 2 -( A +^) t , iX 

H ~ZFt |l)a|l)ba(U| b (U|. 

l + (n) T (l-e- 1 — ) 

I 



Once again we find that at high intensity, i.e., at high temperature, (n)T 3> 1, the last six non diagonal terms 
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FIG. 1: Concurrence as a function of the dimensionless time 
t/tv and of the average photon number of the single mode 
when the initial mode state is: (a) a coherent state and (b) a 
thermal state. In both cases each atom starts in the \0 = 0) 
state. We have considered g/A = 0.01. White color means 
zero entanglement of formation whereas black color stand for 
their maximum entanglement value 1. 



FIG. 2: Concurrence as a function of the dimensionless time 
t/tv and of the average photon number of the single mode 
when the initial mode state is: (a) a coherent state and (b) a 
thermal state. In both cases the two qubits start in a Werner 
state. White color means zero entanglement of formation 
whereas black color stand for their maximum entanglement 
value (37 - l)/2 = 8/11 with 7 = 9/11. 



vanish and reappear suddenly at times t = Like the 
previous case at those tk times the concurrence is zero 
and at any t ^ tk the concurrence is given by Eq. l[8|). 

Figures (p} shows a linear black-white degradation of 
the concurrences of the mixed state given by: (a) the Eq. 
([7]) and (b) the Eq. ([9]), as a function of the dimensionless 
time r = 2g 2 t/ A and of the average photon number of 
the single mode. White color means zero entanglement of 
formation whereas black color stand for they maximum 
entanglement value 1. In both figures ([I]) we considered 
both qubit starting in the \0 — 0) state. 

From the figures ([JJ we see that maximal entangle- 
ment arises periodically at low energy regime. This effect 
is a reminiscence of the entanglement generated by the 
dispersive vacuum state. Those maximal entanglement 
zones are se par ated by narrow entantaglement dead val- 
leys (EDV) [3. 



B. Initial mixed state 

Now we study the case when both atoms are initially 
in a Werner state [2(1 HH type: 

p(Q) = l-ll + 1 \X){Xl (10) 

with I being the identity of the two atoms Hilbert space, 
\X) being one of the four Bell states [2l|, and < 7 < 1 
a physical parameter. 

One can prove easily that when the single mode starts 
in a Fock state \n) the concurrence does not change and 
remains in (37 — l)/2 for 7 > 1/3 and zero otherwise. 

However, when the field starts in a coherent state \a) 
the reduced density operator of the two qubit-system at 



time t becomes 

P = ^/+^(|0) a |0) bQ (0| h <0| + |l) a |l) fcQ (l| h (l| 

±e 2im e-H 2 (i-^)|0) Q |0) bQ (l| b (l| 
±e ^o te -| Q | 2 (i^-^)| 1)a | 1)ba(0 | b(0 | ); (n) 

where we have consider \X) to be one of the two Bell 
states |0±) = (|0>«|0) 6 ± |1)„|1) 6 )/V2. ft = A + ff 2 /A. 
The concurrence of the bipartite mixed state (jlip is given 
by 

fl + 2e- 2 W 2 - 2 ^) 7 -l 
C(P) = ± ^ L , (12) 

for 7 > 1/(1 + 2e- 2 l Q l sin -V) and is zero otherwise. 
Clearly we see that for a high intensity coherent state, in 
general, there is not entanglement between the two atoms 
however at each time t — the initial entanglement 
amount, max{0, (37 — l)/2}, is suddenly recovered and is 
independent of the intensity of the coherent state. When 
one consider the other two Bell states \X) = \ip ) = 
(|0) a |l) b ± |l) a |0) fc )/V2, the (TO]) density operator does 
not evolve. 

On the other hand, when the mode state is initially in 
a thermodynamic equilibrium at absolute temperature T 
the concurrence becomes: 

C(p) = max{0, (^I^Mi-V^)!) 7 " 1 ^ (13) 

Once again we find the entanglement-beats effect at high 
intensity or equivalently in the high temperature regime, 
that is, the initial entanglement amount is suddenly re- 
covered just at each times t = ifc/2. We will call this 
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effect E- beats. We can also seen from Eqs. ([8j and (fl3|) 
that when the field mode is initially in the vacuum state, 
\a\ 2 = (n)T = 0, the entanglement does not change. 
The expression (fTB"]) was calculated considering the state 
\X) = 1^) in Eq. JT[5J). When |A) = the p]) 

density operator does not evolve. 

Figures ((2]) shows a linear black-white degradation of 
the concurrences given by: (a) Eq. (J5J and (b) Eq. (fTB"]) , 
as a function of the dimensionless time r = 2g 2 t/A and 
of the average photon number of the single mode. White 
color means zero entanglement of formation whereas 
black color stand for they maximum entanglement value 
(37 — l)/2 = 8/ll. In both figures ([2]) we considered 
7 = 9/11. From the figures ([2]) we see that the E-beats 
effect becomes apparent even for jMt — \ a \ 2 ~ 3. These 
E-beats are separated by EDVs [151 ] . 

IV. CONCLUSIONS 

In summary we have studied the dynamics of the en- 
tanglement between two non interacting two-level atoms 
weakly coupled and far from the resonance with the same 
single mode field. We find that a dispersive vacuum can 
generate maximum entanglement between them when 
there is a single photon to share. We emphasize that in 
the dispersive regime, the atomic energy is not exchanged 
with the single mode, so the single mode is required to 
be only the mediator between the two two-level atoms 
effective interaction. This effect can not be generated by 
classical field because classical fields can not couple the 
two atoms, at any intensity. When they are initially in a 
type of Werner state, the entanglement is in general zero 



at high energy, but the so called E-beats effect take place 
and the narrow beats are separated by the EDVs [l5j |. 
In other words, in that regime, the initial entanglement 
amount is periodically recoverd in a sudden manner, only 
for a short moments separated by the time scale irA/g 2 . 
The wide of a E-beat is inversely proportional to the en- 
ergy of the single mode. Besides, in that atomic initial 
condition the entanglement does not change when the 
single mode is initially in the vacuum state. However, we 
have already seen that the vacuum initial state makes an 
important effect when each atoms starts in a pure state. 

A physical implementation of this Hamiltonian interac- 
tion between two two-level system and a single mode can 
be performed with two quantum dots interacting with a 
boson mode [22j • Another physical implementation could 
be implemented considing the Zeeman's level structure in 
a 138 Ba cold ion moving in a linear Paul trap [23| in a 
standing wave configuration |24j |. Spontaneous emission 
is suppressed using as a qubit the Sj/2 ground and the 
D5/2 upper metastable states [IH |2(j. The lifetime of 
those metastable states of Ba + is about 45s. The mo- 
tion of the ions can be described in terms of the normal 
center-of-mass mode. The required dispersive interaction 
between two ions with the same center-of-mass mode can 
be always simulated. 
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